Thermodynamics of n-p condensate in asymmetric nuclear matter 
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We study the neutron-proton pairing in nuclear matter as a function of isospin asymmetry at finite 
temperatures and the saturation density using realistic nuclear forces and Brueckner-renormalized 
single particle spectra. Our computation of the thermodynamic quantities shows that while the 
difference of the entropies of the superconducting and normal phases anomalously changes its sign 
as a function of temperature for arbitrary asymmetry, the grand canonical potential does not; the 
superconducting state is found to be stable in the whole temperature-asymmetry plane. The pairing 
gap completely disappears for density-asymmetries exceeding a c = (n n — n p )/n ~ 0.11. 
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Neutron-proton pair correlations are potentially im- 
portant in a number of contexts, including the study of 
the nuclear structure of medium mass N ~ Z nuclei pro- 
duced at the radioactive nuclear beam facilities JlJ and 
the theory of the dcutcron formation in the medium en- 
ergy heavy- ion collisions 0] . In the astrophysical context 
n-p pairingcorrelations are relevant for the astrophysical 
r-process |3|4j and could could play a major role in neu- 
tron star models which permit pion or kaon condensation 

In recent years much theoretical effort has been de- 
voted to the understanding of the n-p pair correlations. 
They have been studied both in the infinite nuclear mat- 
ter within the Thouless or the Bardeen-Cooper- 
Schrieffer || theories of superconductivity |p|-|l6| , and in 
finite nuclei within mean-field effective interaction theo- 
ries of pairing Jl],|l7] . In particular, microscopic cal- 
culations, based on the BCS theory for the bulk nu- 
clear matter show that the isospin-asymmetric matter 
supports Cooper-type pair correlations in the 3 Si- 3 Di 
partial- wave channel due to the tensor component of the 
nuclear force. The energy gap is of the order of 10 MeV 
at the saturation density [B-ll lfl] when the effects of 



the medium polarization on the pairing force 
neglected. 

As the existence of the pair correlations crucially de- 
pends upon the overlap between the neutron and proton 
Fermi surfaces, one expects a suppression of the pairing 
correlations if the system is driven out of the isospin- 
symmetric state. The studies based on the Thouless cri- 
terion for the thermodynamic T-matrix (the diver- 
gence of the of ladder resummation scheme at the critical 
temperature) deduced the suppression of the critical tem- 
perature with the isospin asymmetry fll2| , |l3| |; however, 
they do not permit to draw conclusions about proper- 
ties of the superconducting state with a finite gap. The 
studies based on the BCS theory without self-energy- 
corrections to the single particle spectrum [[l4],[l5|, while 
give a correct qualitative picture, overestimate the mag- 
nitude of the pairing gap and the critical asymmetries at 
which the pairing disappears. 

In this paper we report on a fully microscopic calcula- 
tion of the n-p pairing within the strong coupling theory 



of the superconductivity with realistic nuclear interac- 
tions. Our main focus is the effect of ladder- renormalized 
single particle energies [^9| on the magnitude of the pair- 
ing gap and the threshold asymmetries at which the 
pairing vanishes ]3(| . Our second aim is to provide a 
first computation of the key thermodynamic quantities 
of asymmetric n-p condensate from the strong coupling 
BCS theory. 

To set the stage, let us start with the solution of the 
Dyson equations for the normal and anomalous propaga- 
tors in the Matsubara formalism. In doing so we shall 
decouple isospin singlet SD pairing channel from the 
isospin triple channels, as the SD coupled channels con- 
tain the dominant part of the attractive pairing force. 
In this case the pairing matrix is diagonal in the spin- 
space, i.e. one deals with the unitary triplet state (see 
rcf. uW). The proton/neutron propagators follow from 
the solution of the Gor'kov equations, and can be cast in 
the form (h = 1) 
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8[ are the neutron-proton anomalous propagator has the form 
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where Lu m are the Matsubara frequencies, the upper sign 
in G( p /™) corresponds to protons and the lower to neu- 
trons. The isospin asymmetry lifts the degeneracy of 
the quasiparticle spectra thus leading to two separate 
branches for protons and neutrons 
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and e^ n ' P ' ) are the single particle energies of neutrons and 
protons. The strong coupling BCS theory is coupled to 
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the Brueckner renormalizaton scheme via the single par- 
ticle energies defined as ejj = k 2 /2m + U^ n,p ^(k) — 

^(™iP); here JJ^ n,p \k) are the single particle potentials 
which are derived from the Brueckner theory of asym- 
metric nuclear matter p^] and are the chemical 
potentials for neutrons and protons, which are derived 
from the BCS theory self-consistently. Both schemes are 
normalized to the same densities. The small effect of 
the feedback of the pairing correlation in the Brueckner 
calculations of the mean- field @ is, however, neglected. 

Using the angle-averaging procedure, which is an ad- 
equate approximation for the present purpose (see ref. 
jnj) the BCS gap equation for asymmetric nuclear mat- 
ter can be derived. We find 
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where D(k) 2 = Ao(fc) 2 + A2(fc) 2 is the angle-averaged 
neutron-proton gap function, f(E) = [1 + exp(/3_E)] _1 
is the Fermi distribution function, /3" 1 = fc^T, where T 
is the temperature and ks is the Boltzmann constant. 
The driving term, Vw , is the bare interaction in the SD 
channel. The density matrices of neutrons and protons 
follow from eq. (0]) after summation over frequencies, 
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Summation over frequencies in eq. (0) leads to the den- 
sity matrix of the particles in the condensate 



u(k) = 
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The partial densities of nucleons, in terms of the density 



matrices, are p 



(p/n) _ 



,(p/n) 



(fc); similarly, p° 



Y]^ V(r(k) is the the density of the particles within the 
condensate. It is essential that the coupled system of 
equations (jl]) and (|J) are solved self-consistently. 

Now we are in the position to write down the key ther- 
modynamic quantities. As the occupation of the quasi- 
particle states is given by the Fermi-Dirac distribution 
function, the entropy of the system is given by the com- 
binatorical expression 

S = -2k B k{f(E- k + ) ln/(%+) + f(E n +) ln/(i% + ) 
+ f{E~ k -)\n.f{Ef) + f{Ef) In /(££-)}, (?) 



where f{E^ ± ) = [1 — f(E^)]. Equation (H) is an obvious 
extension of the mean-field expression for the entropy in 
the symmetrical nuclear matter JlT| to the asymmetric 
case. The internal energy, defined as the grand canonical 
statistical average of the Hamiltonian, U = (H— pS"^ p n — 
p^Pp), reads 

U{T,p) = Y^[efn^{k) + efn^\k)_ 

ak 

+ E V{k u -h;k^,-k 2 )v(ki)v(k2), (8) 

where we have carried out the spin summation in the 
second term. The non-pairing interaction energy among 
the quasiparticles and the contribution from the chemical 
potentials are included in the single particle energies (the 
first term) of eq. (||). The second term includes the 
BCS mean field interaction among the particles in the 
condensate. Note that the interaction in the second term 
of eq. (||) can be eliminated in terms of the gap equation 
(H). Finally, the thermodynamic potential is given as 



fl(T,p) = U{T,p)-TS. 
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The expressions for the entropy and the internal energy 
in the normal state follow from eqs. ^ and (||) in the 
limit A = 0. Our main interest below is the change 
in the thermodynamic potential of the superconducting 
state with respect to the normal state 5CI = Q s — fl n . 

Numerical calculations of the pairing gap were carried 
out using the Paris potential. The single particle poten- 
tials were used as an input from Brueckner computation 
based on the Argonne V14 potential; for details of these 
calculations we refer to ref. |2jJ. 

Fig. 1 shows the values of the pairing gap in the 3 Si- 
3 D\ partial wave channel as a function of the temper- 
ature at the saturation density p — 0.17 fm -3 . The 
asymmetry parameter is defined as a = (p n — p p )/p. 
For the value a — one finds the usual BCS solu- 
tion: the gap is a monotonically decreasing function of 
the temperature and vanishes at the critical temperature 
T c = A(T = 0)/1.76. In the asymmetric state the gap 
develops a maximum at a certain intermediate temper- 
ature in the range < T < T c ; the value of the critical 
temperature is reduced and the BCS relation between the 
critical temperature and the value of the gap at T = 
does not hold any longer. For large asymmetries a ~ 0.1 
the superconducting state exists only in a finite temper- 
ature range; it completely vanishes for the critical asym- 
metry a c = 0.11. This behavior is the result of the inter- 
play between the increasing shift between the radii of the 
Fermi spheres of neutrons and protons with increasing 
asymmetry and the smearing of the Fermi surfaces due 
to both interaction driven correlations and the tempera- 
ture. While the first factor tends to suppress the pairing, 
the role of the second factor is twofold. On the one hand 
the temperature and correlations smear the Fermi sur- 
faces thus increasing the overlap between the two Fermi 
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spheres which promotes the pairing. If, however, the tem- 
peratures are high enough the pair correlated states are 
quenched by the thermal excitation. The pairing gap, 
hence, has a maximum at some intermediate tempera- 
ture. The two critical temperatures are controlled by two 
different mechanisms. The superconductivity vanishes 
with decreasing temperature at a lower critical tempera- 
ture when the smearing becomes insufficient to support 
the pairing. The superconductivity vanishes at the upper 
critical temperature because of the thermal excitation of 
the system, as in the standard BCS theory. 

Qualitatively, the role of the mean-field can be under- 
stood in terms of the effective mass. As in the symmet- 
ric matter the nucleon effective mass, m* , is reduced by 
the mean-field, the density of states at the Fermi sur- 
face is reduced, and therefore the pairing gap is reduced 
as well. In asymmetric matter there are no additional 
non-trivial effects. Even though the neutron and proton 
effective masses change with the denisty in opposite di- 
rections, the gap itself depends on the reduced effective 
mass (m* +m* )/2m*m* , which scales with the asymme- 
try as l/m* + 0{a 2 ). 

Fig. 2 displays the entropies of the normal and su- 
perconducting states as a function of the temperature at 
the saturation density. The a — behavior of the nor- 
mal and superconducting entropies is the one expected 
from the theory of the normal Fermi- liquids and the BCS 
theory, respectively: the entropy of the normal state is 
a linear function of the temperature; the entropy of the 
superconducting state is linear close to the critical tem- 
perature and decreases exponentially in the low temper- 
ature limit. For finite asymmetries the entropy of the 
normal state is unchanged. The curve of the entropy of 
superconducting state, however, shows an anomalous be- 
haviour by crossing the entropy curve of the normal state. 
For temperatures below the crossing point the entropy of 
the superconducting state is larger than that of the normal 
state. The crossover to the anomalous behaviour occurs 
at the temperature at which the gap reaches its maxi- 
mum. This anomalous dependence of the entropy of the 
superconducting state on the temperature implies that 
this state could be unstable for temperatures below the 
temperature at which the gap attains its maximum. The 
internal energy, the first term in eq. (^), has a definite 
sign as it vanishes in the normal state and is a quadratic 
form of the pairing gap. Thus the sufficient condition for 
the onset of the instability is that the contribution of the 
second term in eq. ([)]) dominates that of the first in the 
range of the temperatures below the crossover point. 

Fig. 3 shows the difference in the thermodynamic 
potentials of the normal and superconducting state SQ. 
This quantity is negative for all temperatures and asym- 
metries. The anomalous change of sign of S n — S s , there- 
fore, does not change the net balance between the en- 
ergies of the normal and superconducting states. The 
temperature dependence of 8fl is the reminiscent of the 
tempertaure dependence of the gap function by virtue of 



the (approximately) quadratic dependence of the internal 
energy on the pairing gap. We thus conclude that super- 
conducting state is stable, at least in the present model, 
in the whole temperature-asymmetry plane, whenever a 
non-trivial solution to the gap equation exists. We have 
verified that the deduced behaviour of the 5 ft is the con- 
sequence of the fact that the contribution of the internal 
energy to the thermodynamic potential dominates the 
contribution from the entropy. 

The problem of the pairing suppression and the ther- 
modynamics of the superconducting phase when the par- 
ticles lie on different Fermi-surfaces is encountered in 
many systems. The simplest example is a superconduct- 
ing metal in a magnetic field, where the spin degene racy 
is relaxed because of the Pauli paramagnetism [^2|-^4[. 
The spin zero pairing between the spin-up and spin-down 
populations is precisely the one described above, albite 
one commonly uses the BCS weak coupling ansatz which 
does not work for nuclear matter. Another example is 
the B-state of liquid 3 He (Balian-Werthamer state) in 
the magnetic field |3jJ. In both cases measurements of 
the specific heat and spin susceptibility could provide an 
experimental information on the anomalous behaviour of 
the entropy. Apart from condensed matter, our results 
may be relevant for the current discussion of spontaneous 
braking of color symmetry in finite temperature/density 
QCD and formation of (qq) color superconducting con- 
densate (3f] |3Sf1 . In all cases, it would be interesting to 
check whether the balance between the internal energy 
and the entropy is generically such that the supercon- 
ducing state remains stable for arbitrary asymmetries. 

To conclude, we studied the pairing in the isospin 
asymmetric nuclear matter using realistic nuclear interac- 
tions combined with ladder-renormalized single particle 
energies. The pairing at the saturation density vanishes 
for the critical density asymmetry a c = 0.11 at the finite 
temperature T ~ 1.5 MeV. Our evaluation of the ther- 
modynamic quantities of the isospin asymmetric nuclear 
matter shows that, while the entropy of the supercon- 
ducting state becomes larger than that of the normal 
state below a certain temperature, the thermodynami- 
cal potential at the saturation density corresponds to a 
stable superconducting state in the whole asymmetry- 
temperature plane. 
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FIG. 1. The pairing gap at the saturation density as a function 
of the temperature for asymmetries a = 0, 0.05, 0.07 and 0.1. 
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FIG. 2. The entropy of normal (solid line) and superconducting 
state for a = and 0.05, 0.07; the a = 0.1 curve is indistinguishable 
from the normal state entropy curve on figure's scale. 
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FIG. 3. The difference in the thermodynamic potentials of the 
superconducting and normal states at the saturation density as a 
function of the temperature for asymmetries a = 0, 0.05, 0.07 and 
0.1. 
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